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Abstract 

At the LHC, top quark pairs are dominantly produced from gluons, making it difficult to measure 
the top quark forward-backward asymmetry. To improve the asymmetry measurement, we study 
variables that can distinguish between top quarks produced from quarks and those from gluons: the 
invariant mass of the top pair, the rapidity of the top-antitop system in the lab frame, the rapidity 
of the top quark in the top-antitop rest frame, the top quark polarization and the top-antitop 
spin correlation. We combine all the variables in a likelihood discriminant method to separate 
quark-initiated events from gluon-initiated events. We apply our method on models including G- 
prime's and W-prime's motivated by the recent observation of a large top quark forward-backward 
asymmetry at the Tevatron. We have found that the significance of the asymmetry measurement 
can be improved by 10% to 30%. At the same time, the central values of the asymmetry increase 
by 40% to 100%. We have also analytically derived the best spin quantization axes for studying 
top quark polarization as well as spin-correlation for the new physics models. 



1 Introduction 



The recent results on top quark forward-backward asymmetry measurements at the Tevatron have 
shown interesting evidence of new physics beyond the standard model (SM). Large asymmetry was 
observed in both the semi-leptonic decay channel at CDF and DO [2j, and the di-lepton decay 
channel at CDF [3j. While many new physics models have been introduced to explain the top quark 
forward-backward asymmetry [iElElElEllSlinillllliailSlliailllllSllTllig less attention 

has been paid to directly measuring the top quark forward-backward asymmetry ApB at the LHC. 
An obvious obstacle is that the LHC is a proton-proton machine as opposed to the Tevatron, which 
is a proton-anti-proton machine. There is not a universally forward direction for the top quark. This 
obstacle can be overcome by observing that the valence quarks (n or d quarks) are statistically more 
energetic than the sea quarks (n or d quarks). So, event by event one can still define the forward 
direction, making it possible to measure the forward-backward asymmetry [211 122] (see [231 EH ESI [261 
[27] for studies for the LHC). 

Another obstacle is that the tt pairs are dominantly produced from gluon initial states, which serve 
as a huge background for the ApB measurement at the LHC. For the 7 TeV LHC, the production cross 
section of ti from gg in the SM is approximately a factor of 5 (8) larger than from uu (dd). Assuming 
the observed Aps the Tevatron does stem from top quarks produced from up or down quarks, as 
suggested in many beyond-the-SM models, we see that the measurement of Aps at the LHC will be 
diluted by this additional SM background. Although many new physics explanations of Aps will be 
tested indirectly by looking for new resonances in tt or dijet channels, it is still important to have a 
direct measurement of Aps especially when those new resonances are too broad to show up in the 
standard "bump" searches. 

In this paper, we discuss how to improve the Aps measurement at the LHC. More specifically, 
we explore variables that distinguish the uu — t- ti or dd — t- ti production channels from the gg — >• ti 
channel. Because of the differences in parton distribution functions and differential cross sections, the 
produced ti system will have different kinematic distributions. For example, the ti invariant mass 
distribution is often enhanced at large M^i for uu — t- ti or dd — t- ti. The rapidity of the ti system 
in the lab frame, and the top quark rapidity in the ti rest frame also differ for different production 
mechanisms. 

Other than the simple kinematic variables, we also use the top quark polarization and top-antitop 
spin correlation information to reduce the gg — )• ti background. Li the SM, the leading order QCD 
production does not generate polarized top quarks, but this is not the case if the new physics ex- 
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planation of Aps involves new parity-breaking interactions, coupling the left-handed top and the 
right-handed top differently [28l [29l [30l [HI [32] . To study the top quark polarization, one first chooses 
a spin quantization axis and then studies the angles between the top quark decay products and the 
quantization axis. Hence, it is important to know the spin quantization axis. The traditional wisdom 
is to use either the beam direction (beam basis) or the top quark moving direction in the tt center- 
of-mass frame (helicity basis) to quantize the top quark spin. Noticing that neither of those two 
axes can maximize the top polarization effect, we will derive the "best quantization axis" for different 
models and demonstrate its advantages. We emphasize that knowing the "best quantization axis" is 
useful not only for cutting off more SM backgrounds and improve the Aps measurement, but also for 
distinguishing different new physics models. 

Similarly, the spin-correlation between top and anti-top quarks is different for top quark pairs 
produced from gg and from qq. In QCD, when top pairs are produced near threshold, tt is in a ^Si 
state for qq productions and in a ^5*0 state for gg productions [33l[34]- As pointed out in Ref. [35], 
one can calculate the "best quantization axis" to maximize the spin correlation for qq — )• tt in QCD, 
and the so-called "off-diagonal basis" depends on the kinematics of the event. We will follow a 
similar procedure to calculate the "best axis" in the presence of new physics contribution to qq — )• ti, 
especially for those models without polarized tops. The top quark spin-correlation is an important 
effect to distinguish the tt resonances as studied in Ref. [361 [33 [Ml [39] ■ We believe that the formulas 
developed in this paper will be useful for identifying the tt resonance properties once the LHC has 
positive results in those searches. 

We combine the variables in a multivariate likelihood method to improve the significance of the 
asymmetry measurement, defined as Afb/(^Apb- We find that the significance can be improved by 
10% to 30% for the models we consider. Meanwhile, the central values of the asymmetry increase by 
40% to 100%, making the measurement less sensitive to systematic uncertainties. 

Our paper is organized as follows. In section [2l we choose three representative models that explain 
the ApB results at the Tevatron and calculate their predictions oi ApB for the LHC without using 
additional variables. In section [31 we consider three basic kinematic variables: the ti invariant mass. 
Mjj, the boost of the system, y^j, and the rapidity of the top quark in the rest frame, yt- We study 
differences between gg — t- tt and qq — t- tt in terms of those variables. We then study top polarization 
in section [D and top-antitop spin-correlation in section \5\ with more details given in appendix [3 In 
section [6l we demonstrate the combined improvements using all the variables. We then discuss various 
effects such as experimental cuts and event reconstructions on the improvements when performing a 
realistic analysis, and conclude our paper in section [3 
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2 Top quark forward-backward asymmetry at the LHC 



To explain the large top quark forward-backward asymmetry measured at the Tevatron, there are 
two basic top quark pair production mechanisms. One is through new particle exchange in the s- 
channel such as the axigluon, G' , and the other one is through new particles in the t-channel such as 
the diquark, Z' or W'. Noticing that Z' model generically predicts copious same-sign tops, which is 
tightly constrained by the recent analysis from CMS [40], we choose G' and W as two representative 
examples in this paper. We will also consider effective contact operators obtained by integrating out 
a very heavy axigluon. G' and W couple to the SM quarks as 

Cg' = -G'; [n (g^7^i" + 9^'^'^^ + ^(5^7^*" + 5^7^'*")*] + • • • , (1) 
= -W'+ i{gvl^ + 9Alh'')d + h.c. + • • • , (2) 

where we only write down couplings relevant to Aps and is the SU{3)c generator. The differential 
production cross section as a function of the top quark production angle is given in [41^ |3] for the 
axigluon case and in [42] for the W case. The measured Aps Si-t the parton level is [l] 

ApBiMtt < 450 GeV) = -0.116 ± 0.153 , ApBiMu > 450 GeV) = 0.475 ± 0.114 , (3) 
^Fij(|Ay| < 1.0) = 0.026 ±0.118, ^Fij(|Ay| > 1.0) = 0.611 ± 0.256 . (4) 

To fit those experimental data, we consider the following three different model points: 

• Model A: an axigluon model with Mc = 2.0 TeV, g\ = 2.2, 5^ = -3.2, g{, = 1.0 and = 0. 
Here, "g" represents the first four light quarks. The width is Tc' = as/6{4g'^ + Agy^ + 2g^^ + 
2g^) Mc ~ 1.5 TeV. The predictions for ApB are (0.10,0.31) for the two invariant mass bins 
and (0.12,0.40) for the two rapidity difference bins. 

• Model B: a W model with Mw' = 400 GeV and gy = gA = 0.9 (or gi = and gR = 1.8). The 
predictions for Aps ai^e (0.12,0.41) for the two invariant mass bins and (0.14,0.52) for the two 
rapidity difference bins. 

• Model C: the contact interaction obtained by integrating out a very heavy axigluon (above the 
center of mass energy of the LHC and the Tevatron), S, u'j^'y^t'^ut^'^^^t^'t / A'^ , with g^* = 0, 
^ = -1 and A = Mq' /{g^g'x)^/'^ = 650 GeV. The predictions for ApB are (0.19,0.53) and 
(0.23, 0.67) for those four bin data. 

All of the above models satisfy various constraints such as those from dijet resonance, dijet contact 
interaction and ti resonance searches at Tevatron [9] , though they receive more stringent constraints 
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from the latest LHC results on tt resonance searches [43] . Model A has an axigluon coupling to the 
top quark with both vector and axial-vector couplings, so the parity is broken. As we will show later, 
the top quark will be polarized in this case and we can use the top quark polarization to reduce the gg 
background and improve the Aps measurement. The axigluon width is very large for this model, so 
the tt resonance would be difficult to discover. For Model B, we only choose right handed couplings 
for the W' field, based on the severe constraints from electroweak precision observables if W mixes 
with the SM W gauge bosons. The W width is assumed to be small for this case and neglected in our 
later analysis. Model C only has parity conserving couplings and will not have polarized top quarks. 
We use this model as an example to study the top and anti-top spin-correlation. 

The measurement of top quark ApB at the LHC is challenging for two reasons. First, it is a 
proton-proton collider. Unlike Tevatron with proton-antiproton collisions, there is not a fixed forward 
direction. However, the valence quarks inside the proton most likely carry a larger energy than the sea 
quarks. So, event by event, we can define the moving direction of the center-of-mass frame with respect 
to the lab frame as the positive direction to calculate Aps- There is a ~ 20% probability that we 
misidentify the initial parton moving directions. The second reason is that the main mechanism for tt 
production is gg — )• tt, which does not contribute to Aps- Using Madgraph we obtain the leading 
order QCD cross sections as cr{gg — )• tt) « 71 pb, a{uu — t- tt) ~ 14.5 pb and a{dd — t- tt) ~ 8.6 pb for 
the 7TeV LHC. One can immediately see that there is an additional tt background from gg by a factor 
of 5. This is different from the Tevatron case, where cr{gg tt) is a subdominant part in the total 
production cross section. Including all tt pairs and neglecting other SM backgrounds, we follow the 
above definition of the forward direction and estimate the Aps measurements for the three models as 

Model A: ApsiMn > 450 GeV) = 0.046 ± 0.015 , 

Model B: ApsiMn > 450 GeV) = 0.196 ± 0.011 , 

Model C: Apb{Mu > 450 GeV) = 0.099 ± 0.015 , (5) 

where, to estimate the statistic errors we have assumed a 3 fb~^ luminosity and semileptonic decays for 

the tt system. Furthermore, we have multiplied the total number of events by a 10% event acceptance 

(depending on cuts, the acceptance could be even higher than this value, see Ref. [l3] for example), 

so there are around 5000 events in total. From the numbers in Eq. ([5]), one can already see that the 

early LHC running can measure Aps at a large confidence levej^. 

^After imposing the Mu > 450 GeV cut, the production cross section for gg — >■ ti is 37.6 pb, the "signal" production 
cross sections are 13.1, 46.6 and 17.9 pb for Model A, B and C, respectively. Model B predicts too many tt's and 
has already been ruled out by the M-tt differential cross section distribution [43] • Here, we still include this model to 
illustrate how to improve the LHC Afb measurement for a t-channel model and we will also consider a similar model 
with a smaller coupling. 
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To improve the measurement of Aps at the LHC, we need to distinguish the two tt production 
mechanisms: from two gluons or from two hght quarks. In the following sections, we will study the 
differences between various kinematics distributions of tt for those two production mechanisms. 

3 Basic kinematics: invariant mass, M^f, the rapidity of the ti sys- 
tem, y^f, and the rapidity of the top quark, yt 

In this section, we first consider some basic kinematic distributions for events generated from the 
process gg — )■ tt and from the three model points. 

The first variable we consider is the tt invariant mass. Since there are new heavy particles con- 
tributing to the top pair production, we anticipate that the tail of the differential cross section in terms 
of should be lifted. Therefore, the top pairs from new physics should have a harder spectrum than 
from gg. This is illustrated in Fig. [1] (left panel). Note that Model B has a much harder spectrum than 
the other cases. So, we can impose a higher M^j cut to enhance the signal-background ratio. However, 
as we will show in Section [6l cutting on Mjf alone may not increase the statistical significance and 
may not improve the Aps measurement. Therefore, M^j should be combined with other variables. 




Figure 1: Left panel: the normalized fraction of events as a function of M^f. Model A and Model C 
contains both uu — )■ tt and dd — )• tt productions, while Model B only contains dd — t- tt. Right panel: 
the rapidity distributions of the center-of-mass frame of the ti system. 



The second variable is the boost of the ti system with respect to the lab frame y^f. From the 
Parton Distribution Functions (PDF's), one expects from uu productions to be statistically larger 
than from gg. This is indeed the case as shown in the right panel of Fig. [1] from which one can see 
that Model A and Model C have more signal events at larger values of \ytt\- 
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Next, we consider the rapidity of the top quark in the tt center-of-mass frame: \yt\- This variable 
is especially useful for selecting signal events for Model B. This is because of the t-channel differential 
cross section enhancement in the forward direction. From the simulated result in Fig. [2l we see that the 
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Figure 2: The rapidity distributions of top quark in the center-of-mass frame of the tt system. 

\yt\ distribution in Model B (the dot-dashed blue histogram) peaks at around 0.8, which is significantly 
different from the background (the solid black histogram). On the contrary, for Model A and Model 
C, the produced top quarks are a little more central but similar to the background. So, we do not 
anticipate \yt\ as a good variable to improve the Aps measurement for Model A and Model C. 

4 Top quark polarization 

The top quark has a short lifetime and decays before hadronization, so its spin information is kept 
in the angular distributions of the daughter particles. Its polarization as well as spin-correlation with 
the anti-top quark are different for new physics and the SM. In this section, we utilize the top quark 
polarization to distinguish signals from the gg background . 

Due to parity conservation, top quarks produced from QCD processes are not polarized. For Model 
A and Model B, parity is manifestly broken and top quarks generated from new physics are polarized. 
Choosing a spin-quantization axis for the top quark, one can study the angular distribution of the 
daughter particles in the top quark rest frame, with respect to the axis. The differential decay rates 
for a 100% polarized top in its rest frame are calculated in Ref. |28j 






(6) 
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where 9i is the angle between the chosen spin-quantization axis and the i'th decay product in the top 
rest frame. For top quark, one has k£+ = = kg = I, ky^ = ku = k^ = —0.31, k\) = —k-^+ = —0.41. 

Now, the important task is to find the best spin-quantization basis to keep the maximal top 
polarization information. There are two obvious basis: the beam line basis (the moving direction of 
the incoming quark in the top rest frame) or the helicity basis (the top quark's moving direction in 
the ti rest frame). The spin density matrix can be decomposed into independent parts with different 
spin structures: 

p = yll2 <^ 12 + Bt • o- (g) I2 + Bf • I2 ® o- Cijo' (g) , (7) 

where Bj (Bj) determines the top (anti-top) quark polarization and Cij denotes the potential spin- 
correlation between top and anti-top, which is especially useful when = Bj = 0. The detailed 
calculations and formulas for those vectors and matrices can be found in Appendix lAl The polarization 
vector B( can be decomposed into 

Bi = 6fp + 6fk. (8) 

Here, p = (0,0, ±1)-^ is the beam direction and k = (0, sg*, c^*)^ is the direction of the top quark in 
the ti rest frame where 9* is the top production angle. We have assumed CPT is conserved for the 
underlying theory, so the projection on the p x k direction vanishes. Furthermore, we also assume CP 
is a good symmetry such that 6^ = 6^ and = b^. Therefore, the anti-top has the same polarization 
as the top, and we only need to consider the top quark polarization in this section. 

Starting from the axigluon case. Model A, with 5^ = and gy 7^ 0, we use the following ratio to 
define the best quantization basis for top quark polarization: 



b 



k Acose*{Mli,-s){2mt-\fl)+g\g\sj3 cos2e*{\fl - 2mt) - 2mt + 



bP Amt{g\g\sP cos e* -2M I,, + 2s) 



(9) 



where /3 = y^l — Amf/s. One can see that when /3 = or \/I = 2mt, the above ratio is zero. This 
means when two top quarks are produced at rest, the beam line is the best quantization basis. In the 
other limit when s = Mq, (or when the axigluon is on-shell), this ratio is (sec 6* + cos 9*)Mg' / {2mt) — 
cos 6*. So, for Mqi ^ nit the best spin quantization basis is the helicity basis k. We also note that 
when gy = gy = the top quark is not polarized. In the left panel of Fig. [3l we show ratios of b^ 
over 6j as a function of s for different production angles. For s between 500 GeV to 1 TeV, 6j and 
are comparable to each other, then the best spin quantization axis is neither the helicity basis nor the 
beam axis. So, using the best spin quantization axis may significantly increase the top polarization 
measurement and help improve the ApB measurement. In the right panel of Fig. [3l we compare the 
sizes of the top quark polarization for different spin quantization axes. One can see that the top 
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Figure 3: Left panel: the ratio h^/lfl. of Model A, which determines the best top spin quantization 
axis, for different kinematics. Right panel: the normalized distributions of the angle between the 
lepton momentum in the top rest frame and different top spin quantization axes. All momenta are 
first boosted into the ti zero momentum frame and only the signal events, ti produced from light 
quarks, are included in this plot. 



quark is indeed polarized in Model A. Since the background top quarks from gg are not polarized, 
the corresponding cos0£+ distribution is flat. Therefore one can use polarization effects to distinguish 
signals from backgrounds. From this plot one also sees that the "best quantization axis" defined in 
Eq. ([9]) gives us a larger polarization effect than the other two axes: the helicity basis and the beam 
line basis. The differences could be even larger assuming we know exactly the u or d parton directions 
in each event. Unfortunately, the LHC is a proton-proton collider and there is a ~ 20% probability 
that we misidentify the initial parton moving directions. 

Let us turn to Model B. The W boson in Model B only has right-handed couplings to quarks. A 
large polarization effect is anticipated in this model. Similar to the Model A case, we use 5j and 
to define the best quantization basis for top quark polarization: 

- 8m'^M^,s [9(3glmsz + 2t (2m + {(3z - if)] - 32M^, [2m + {f3z + if] 

+4M^, [9gjis{Pz + l)(2m + f3z + l) + 8m^s + im'^s{f3z - if + 16mt + 8t{pz + if] } , 

2 

bf = ^ X {-9glm'^s^((3z-l) + 4.m'^M^,s(9gls-4:l3tz + 8t) 

-32M^,{(3z + 2) + 4M^, (9g|(/3sz + s) + 4m'^s{/3z - 2) + 8f3tz + IQt) } . (10) 
Here, t = —js{l + — 2/3z), z = cos 9* and m = mt/\/s. The right-handed gauge coupling gji is 
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normalized with respect to the QCD couphng gs- To derive those formulas, we have neglected the W 
widths. The ratios /b^ as a function of the s for different production angles are shown in the left 
panel of Fig. HI The comparison of the top polarization for three different spin-quantization axis is 
shown in the right panel of Fig. [H One can see that there is not much difference between the "best 
axis" and the helicity basis. The reason is that in this model most top pairs are produced with large 
center-of-mass energies, which can be seen from Fig. [TJ Because of the chiral coupling, in the massless 
limit the top quark has a definite helicity, i.e., it is 100% polarized in the helicity basis. Therefore, 
the large momentum of the top quark makes its mass unimportant and the helicity basis close to the 
"best axis". The top quarks are largely polarized for Model B as can be seen from the right panel of 




(GeV) COS0,* 

Figure 4: The same as Fig. [3l but for Model B. 



Fig. m Therefore, imposing a cut on cos6'£+ may improve the ApB measurement a lot. 

5 Top and anti-top quark spin correlation 

Top quarks are not polarized in Model C, but the spins of the top and the anti-top are correlated. In 
this section, we show how one can obtain the correlation information and use it to distinguish between 
events from gg and from qq. 

The spin correlation tt can be measured by studying the following double differential distribu- 
tions 

1 (fN 1 

-{1 — C kikj cos 9i cos 9j) , (11) 



N d cos 9i d cos 9j 



The angles 9i {9j) is the angle between the quantization axis in Eq. ()13p and the daughter particle's 
momentum from top (anti-top) decay, measured in the top (anti-top) rest framed. The coefficients ki, 
^One should first boost every momentum into the tt rest frame, and then boost the particles in the final state into 
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kj are constants determined by the particle species [33]. Instead of fitting the distribution in terms 
of two variables cos 9i and cos 9j , one can integrate Eq. ()lip to obtain the following one-dimensional 
distribution 

1 dN 1 „ . 

(12) 



dfcos 9-1 cos ( 



- (C kikj cos 9i cos 9j — 1) log(| cos 9i cos 9j \ 



The parameter C depends on the spin-quantization axis. The "best axis" to maximize the spin 
correlation for Model C is (for the detailed derivations, see Appendix I A. 2 p 



oc p + 



ce*{-y - 1) 



A2 



k, 



(13) 



In the left panel of Fig. \5\ we show the ratios of the k component over the p component for different 
s and 9* . In the right panel of Fig. EJ we compare the spin-correlation effects by using different spin- 
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Figure 5: Left panel: the same as the left panel in Fig. [3l but for Model C. Right panel: the normalized 
event distributions in cos^£+ cos^^- for different spin quantization axes. 



quantization axes. One can indeed see that using the "best axis" can increase the spin correlation 
compared to the other two axes. If there is no spin-correlation, the distribution should be symmetric for 
positive and negative values. Using the analytic formula in Eq. (|12p . we fit the simulated distributions 
and found that C = 0.47, 0.80, 0.94 for the beam line, helicity and "best" axes, respectively. 

The background events from gg — )■ ti should also have spin-correlations for those three spin- 
quantization axes. In Fig. [6l we compare the signal .vs. background distributions using the "best 
axis". One can see that their distributions are indeed different from each other, which can be used 
to improve the ApB measurement. For the di-lepton channel, we can use the two charged leptons to 
study the spin-correlation. After fitting the distributions, we have C = —0.13 for gg — t- ti, which has an 
the top/anti-top quark rest frame. 
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Figure 6: Left panel: the normalized event distributions of cos0£+cos^£- for Model C and the tt 
events produced from gg in the di-lepton channel by using the "best axis" . Right panel: the same as 
the left panel but for the semi-lepton channel. 



opposite sign to the signal, C = 0.94. For the semi-leptonic channel, we can identify the jet closer to the 
6-quark in the gauge boson rest frame as the down-type quark (the probability is around 60% from 
Ref. [33]). After fitting to the distributions, we have C = —0.08 for gg — )• tt, which has an opposite sign 
to the signal, C = 0.41. Since gg — t- ti dominates the tt productions at the LHC, the spin correlation 
from the signal model will be diluted. making it difficult to measure the spin-correlation. However, 
eventually the LHC may accumulate enough data and make the spin-correlation measurement feasible, 
then one needs to find the optimal spin-quantization axis for the gg — t- tt productions and the formulas 
obtained in this paper would be useful not only for the measurement but also for distinguishing between 
models. 

6 Combined Improvement 

In this section, we combine the useful variables defined in the previous sections and consider the 
improvement on the Aps measurement. Instead of using simple rectangular cuts, we adopt a likelihood 
discriminant method described as follows [45]. 

For a given variable x*, we obtain from simulation the signal and background distributions as given 
in histograms s* and 6*. We normalize and 6* such that they have the same binning and area. For 
a given event with the variable falling in the j'th bin, we define the probability of it being a signal 
event as 
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where s*- and 6*- are the numbers of events in the j'th bin for histograms and 6* respectively. For 
multiple variables, the signal likelihood is defined as 

Next we need to specify what variables to use when calculating the likelihood. The kinematic 
variables, M^j, \yt\ and are useful for all models so we will always include them. Model B has 
large polarizations for the tops, so we include cos 6e and cos 6d using the best quantization axis for 
Model B, but not the spin correlation variables. The tops are not polarized in Model C, so we will 
use the spin correlation variable cos 0^ cos 0^ based on the best quantization basis. For Model A, we 
use both the polarization variables and the correlation variable. For all variables, we group them to 
histograms with 20 bins. For cosOi, cos 9 a and cos 9i cos 9^, the 20 bins have the same size from -1 to 
1. For M^f, the first 19 bins have a size of 25 GeV ranging from 450 GeV to 925 GeV and the last bin 
contains all events with Mjf > 925 GeV. For \yt\ (|yttl); the first 19 bins correspond to (0, 1.9), evenly 
distributed, and the last bin contains all events with \yt\ > 1.9 {\ytt\ > 1.9). 




Figure 7: Likelihood distributions for signal and background events. 



The signal and background likelihood distributions for the three models are shown in Fig. [71 
Given the distributions in Fig. [TJ we can choose a particular likelihood cut Cf^^ and keep events 
with Cs > Cf^^. We then obtain the significance, S/\/S + B, as a function of >C™*- Since the signal 
efficiency is a monotonous function of the likelihood cut, we can change the variable and draw the 
the significance as a function of the signal efficiency instead, which is shown in Fig. [HI Note that the 
peaks at high likelihood values arise from the high Mjf tail, where the events are much more likely 
signal events. The most efficient variable to distinguish signal and background is also M^j. Therefore, 
it is illuminating to examine the significance improvement by cutting on M^f alone and to compare it 
to the improvement by including all variables, which is also shown in Fig. [HI Another characteristic 
of the improvement is the signal-background ratio (S/B). Obviously, for the same significance, we 
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Figure 8: Significance {S/ \/ S + B) as a function of the signal efficiency. 

would like S/B as large as possible. This is due to two reasons. First, for larger S/B, the results will 
be less sensitive to the background systematic uncertainties. Second, since the background in our case 
has no asymmetry, for larger S/B, the measured central value of asymmetry will also be larger and 
deviate more from a flat distribution. From Fig. [8] and Fig. [9l we see that we can increase S/B and 
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Figure 9: Signal-background ratio as a function of the signal efficiency. 

at the same time obtain moderate improvement in S/\/ S + B. This will help with the discovery of 
the signal events by making it less sensitive to systematic uncertainties. 

On the other hand, we should make sure that the cuts we use do not significantly reduce the 
asymmetry. Therefore, we need to calculate the asymmetry using events which pass the likelihood 
cut and examine directly whether we have improved the asymmetry measurement or not. Moreover, 
the likelihood cut that maximizes S/\/S + B in general is not the cut that maximizes the asymmetry. 
Therefore we scan the likelihood cut and find the cut that maximizes Afb/(^Apb^ where ctafs is the 
error for the asymmetry measurement given by 1 /^fN with the total number of events after cuts 
and after taking a 10% acceptance into account. The asymmetries and errors, as well as their ratios 
as a function of signal efficiency are shown in Fig. [101 Then we find the best likelihood cuts for the 
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Figure 10: Asymmetry central values and errors as a function of signal efficiency. 



three models are 0.49, 0.33 and 0.36 respectively, corresponding to signal efficiencies of 0.60, 0.80 and 
0.86, and background efficiencies of 0.30, 0.37 and 0.63. The resulting asymmetries are given by 

Model A: Afb{Cs > 0.49) = 0.078 ± 0.024 , 

Model B: Afb{Cs > 0.33) = 0.289 ± 0.014 , 

Model C: ^fbC-Cs > 0.36) = 0.121 ± 0.017 . (16) 

Comparing with the numbers in Eq. ([5]), we see that we have achieved larger central values for the Aps 
measurements with improved ApB /(^ApB ■ Note that the improvement in Aps/ (^Aps significant 
for all three models. This is due to different reasons: for Model A and Model C, the likelihood 
distributions for the signal and the background are not dramatically different; for Model B, although 
the distinction between the signal and the background is large, the signal cross section for m^j is so 
large that it is not essential to reduce the number of background events. As mentioned previously, 
given the recent LHC results on Mjf distribution measurement. Model B is no longer viable unless the 
W' — t — d coupling, and therefore the signal cross section are smaller. In that case, our method will be 
more useful. As an illustration, we consider Model B with the same W mass but a coupling gpt^ = 1.5 
and repeat our optimization procedure. The signal cross section for M^f > 450 GeV is reduced to 21 
pb from 38 pb of the original model. The resulting asymmetry as a function of signal efficiency is 
given in Fig. [TTl We see that ApB / oapb is improved by about 30% for the best cut, with the central 
ApB value more than doubled. 
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Figure 11: Asymmetry central values and errors as a function of signal efficiency for Model B with a 
smaller coupling, gn = 1.5. 

7 Discussion and conclusions 

For all new physics explanations oi ApB, the differential cross sections in have been predicted to 
be different from the SM. Measuring this distribution would be the first hint of new physics behind 
ApB- However, the measurement of the production cross section as well as da/dM^j; distributions 
are suffering from a large systematic errors related to the jet energy scaling and the luminosity uncer- 
tainty [43|. If the resonance in the s-channel is very broad or the new particle in the t-channel does 
not contribute to the tt productions significantly, performing a precision measurement like measuring 
ApB could be the unique way to unravel the new physics behind top quarks. 

Treating the tt production from gluons as backgrounds and those from light quarks as signals, we 
have found that the central values oi Aps can be increased by a factor as large as 2 for all three 
models when we impose stringent cuts on both signal and background events. The real improvement 
on the significance of measurement, on the contrary, can only be increased by 10% to 30%. The simple 
reason is that the optimized cuts from our likelihood analysis decrease the signal efficiency as quickly 
as increase the central values Aps- We believe that this result is true not only for the measurement 
of ApB but also for the charge asymmetry measurement because of a strong correlation between ApB 
and the top quark charge asymmetry. From Fig. [9l one can see that the ratio of top quarks produced 
from light quarks over from gluons can be increased by 100% with a moderate cost of signal efficiency. 
This improvement of S/B can eventually help the ApB measurement especially when the systematical 
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errors are large. 

In all of our analysis above, we have neglected the effects on our results from the top reconstruc- 
tion procedure as well as various experimental cuts on the final state particles from the top decays. 
Especially, the top polarization and top spin-correlations will be affected by the rapidity cuts on lep- 
tons and jets. The actual difference between the signal and background distributions in Fig. [6] will 
be reduced. To find more accurate experimental improvement in the Afb measurement using the 
variables in this paper, one also needs to consider other effects such as hadronization, initial state 
radiation, final state radiation and the detector effects, which are not included in this paper. 

Another concentration of this paper is to find the best quantization axis for studying the top quark 
polarization and top-antitop spin-correlation. For the new physics models considered in this article, 
we have obtained the axis giving us the maximal top polarization, as well as the axis yielding maximal 
top-antitop spin correlation. If new particle contributions to top quark pair production do exist, one 
can scan different spin quantization axes to find the best one that maximizes the polarization or spin- 
correlation effects. Then, we can use the analytical formulas for the best quantization axis to fix a 
relation among the model parameters. From this point of view, finding those axes are not only helpful 
for improving ApB measurement but also useful for disentangling different new physics models. 

In conclusion, we have explored various kinematic distributions for improving the Aps measure- 
ment at the LHC. We have adopted a multivariate likelihood discriminant to obtain the potential 
improvement by including the invariant mass of the top pair, the rapidity of tt system in the lab 
frame, the rapidity of top quarks in the ti rest frame, as well as top polarization and top-antitop spin 
correlation. Treating top pairs produced from gluons as the background and those from light quarks 
as the signal, we have found that the Aps measurement can be improved by 10% — 30% depending on 
the underlying models. The ratio of tt production from light quarks over from gluons can be increased 
by as large as 100%. We have also included our calculated best spin quantization axis to maximize 
the top quark polarization and spin-correlation effects in the appendix. We believe that those axes 
could be very useful to distinguish different new physics models eventually. 
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A Formulas for top polarization and spin correlation 

A.l The axigluon case 

We assume that the axigluon has universal vector and axial couplings to the first two generations, gy 
and g'^, while the couplings to the top quark and bottom quark will be gy and (7^. For convenience 
we have rescaled these couplings in terms of the strong coupling gg. The tree- level differential cross 
section for qq — )• tt will be |4H S] 



dcosO* ~ 9J 



, 4 2 , 2^ 2glg\,s{Ml,, - s) gVjgV + g\^ 

^ M (S - M2,)2 + M2,r2. ^ (s - M2,)2 + M2,r2. 



S2 



„ J / s(M^, - s) q ^ 52 

" ^^^^'^ I^(s-M2,)2 + M2,r2, " ^^^^^(s-M2,)2 + M2,r2, 



(17) 



where = ml/s, /3 = \/l — 4m2, c = /3cos0*, and 0* is the angle between the top quark and the 
incoming quark in the center of mass frame. The forward-backward asymmetry arises solely from the 
last line in this equation, so to obtain a positive asymmetry we must have g\g\ ^ ^• 

To study the top and anti-top quark polarizations and the spin correlations of top and anti-top 
quarks, we calculate the spin-density of the top quark pair productions, following Ref. [l6] and Ref. [37] . 
The matrix element square is decomposed into different spin structures of top and anti-top 



22iV2 

spin{i'nitial) , color 



(18) 



Here, st (sf) is the unit polarization (three-component) of the top (anti-top) quark in the rest frame 
of the top (anti-top) quark and cr* are Pauli matrices. The spin density function p is a 4 x 4 matrix 
in terms of the four-components of the spins of top and anti-top, which defined as 

k • St k(k • Si) 

St + 



mt mt{mt + E) 

= ';"'•:'' ). (19) 

V "^i rntimt + E) J 

Here, E and k are the energy and the three-momentum of the top quark in the ti rest frame. Those two 
four-component spin vectors satisfy the Bjorken and Drell relations: sf = s^ = — 1 and Sfkt = Sf-kf = 
where and A:^ are the top and anti-top quarks four-momenta in the center-of-mass frame. To 
calculate the matrix element, the following relations of spinor outer products uu and vv are useful: 
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u{p, s)u{p, s) = ^(2/ + m)(l + 75^ and v{p, s)v{p, s) = ^(^Z — m)(l + 75^) (see Ref. [33] for more detailed 
discussions). 

The spin density matrix is decomposed into independent parts with different spin structures 

/) = AI2 ® I2 + Bt • (T » I2 + Bf • I2 ® + Cijd' (g) , (20) 

whether B^ (Bj) determines the top (anti-top) quark polarizations and Cij denotes the potential 
spin-correlations between top and anti-top, especially when B^ = Bj = 0. To determine the spin- 
quantization axis and without loss of generalities, we choose a coordinate with the unit vectors of the 
initial up quark and the final top quark as 

p = (0,0,1)^, k={0,sg,,ce*f , (21) 

with So* = sin 9* and cg* = cos 9*. One can also define one more unit vector perpendicular to both p 
and k as n = p X k/|p x k|. The matrices Bj, Bf, C can be further decomposed into 

Bt = hfp + blk + bfh, (22) 
Bf = 5fp + fefk + 6fn, (23) 
C = col3 + C4P0p + C5k(g)k-Fc6(k0p-Fp(g)k). (24) 

Here, assuming CPT is a good symmetry, we have = b'^ = 0. Furthermore, we will assume that CP 
is a good symmetry for the axigluon model, so one has 6f = 6| and b^ = b^. Other structures of the 
matrix C are also forbidden by CP and CPT symmetries. Under P-symmetry transformation, one has 
ti^^ — )• — ^£ and 6^. — )• — So, if parity-symmetry is a good symmetry, we anticipate B^ = Bf = 0. 
In the QCD and for the production processes qq — t- tt at tree level, one has 

^QCD ^ L^2-I3'sl), B,Q^° = bQ^° = 0, (25) 
18 

QCD _ ■Sg* QCD _ 1 

•^0 — "l^g ' "^4 — g ' l-^OJ 

(9m — ^^r■^^ 

(27) 



QCD _ (4m - phi + 2) QCD _ (2m - 1) cg 



5 9(2m + l)2 ' b 9 • 

An overall coefficient g'^ is understood. The interference terms of QCD and axigluon contributions 
have the following formulas 

s (s — ) 

= i8[(i-Mg,)^°Mg,ra.] [^•"»^»^^«- + ^•"'^''a ■ 

S (s — ) 

^ 18[(s-M2,)2 + M2,r|,] L ^ Ayi/yAM » yyyA;j > v ; 
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and 

[-2g^g\rP\sl] , (31) 

[^aWv] , (32) 
[Ag\g\{2m + 1)/?^ ce> + 4/3^ gt-ff|^(4m - P^s^ + 2)] 



Jnt _ 




s (s — 


Ml,) 


— 


18 [(s 




)' + Ml,rl,] 


Jnt _ 
C4 — 




s [s — 


Ml,) 


18 [(.§ 




)' + Ml,Tl,] 


jnt _ 




s [s — 


Ml,) 


— 


18 [(s 


-Ml,: 


)' + Ml,rl,] 


Jnt _ 




s [s — 


Ml,) 



(33) 



(2m + 1)2 

•(2m-l)ce*l, (34) 



18[(s-M2,)2 + M2,r2.,J L ™ . ■ .v.v. 

The new-physics-only part has 
.new ^ {\{9^^ + ffl^) [/3'(gt-^ + 9f)c2e* + 3/3^^^^ + gt-^(4 - p^)] + Sgl^fft^g^^g^/? cg. } 

18[(5-M2,)2 + M2,r2.,] ' ^""^^ 

,new ^ [4 m fft. g\ PiaA" + + 8 ^ 9\9Uv^\ 

P 18[is-Ml,)^ + Ml,Tl,] ' ^ ^ 

- 2 /? §2 

18(2m + l)2s0* [(s-M2,)2 + M2,r2,,] ^^^^^^ V J yy y Jl 

+ dv'W [(2m + l)/3'c2e. + (2m - l)/?^ + 4(2m + 1)] | . (37) 



k 



and 



18 [(S-M2,)2+M2,r2,,] 



(38) 



new ^ [2{gl- + gl-){g\^ - g\^)] 

' 18[(5-M2,)2+M2,r2,,] ' ^""^^ 



5' [SI5h\g\,g\g^{2m + 1)0^3. + /^^^^'(^l' + 5^>')(8m - 2/32^2. + 4)] 



18 [(s - M2,)2 + M2,r|,] (2m + 1) 



2 



(40) 



new ^ [2{gf +gX.''W^P^ + g\^{2m-l))ce, +Sg\glg\,g^m0\ 

18 [{s-Ml,y + Ml,Tl,] ' ^ ' 

From Eq. (jl7p and to maximize the effects of increasing Ap^, we choose 9y = 5y = 0. For this 
choice of parameters, ah 6's are zero and there is no polarizations of top and anti-top quarks at the 
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tree level. The summations of c's are 

o2 „2 -2o2 t2 u2 2 

" 18 ^ 18 [{s - M2,)2 + M2,r^,] ' 

1 -2o2 u2 t2 

' 9 9[(s-M2,)2 + M2,r|,] ' ^ 

(4m - /3^,s^, +2) s {s - M^,) 2glg\(3^ cp* 

' 9(2rn + l)2 ^ 9 [(s - M2,)2 + M2,r^,] (2m + 1) ' ^ ^ 

(2m-l)cg. g(g-M^,)2ff^g>/3 /g^'gAV 

9 9[{s- M2,)2 + M2,r2,,] ^ 9 [(s - M2,)2 + M2,r2,,] ' ^ 

To determine the spin axis to maximize the spin correlations between top and anti-top quarks, one 
needs to diagonalize the Cij matrix 

CO 

C= I co + c^sj, CQSg* + C5Ce*sg* | , (46) 

CQSg* + C5Cg*Sg* Cq + C4 + C^Cg, + 2cqC0* 
and find the eigenvector corresponding to the largest eigenvalue. 

A. 2 Contact Operator 

We first analyze the contact operator case. The relevant operator can be written as 1x7^ r" j^u tj'^j5 r" i/A^ 
with SU (3)qcd generators. To match the formulas for the axigluon model, we identify 



M2, A2 ' 



(47) 



and Tg' = 0. Here, ^ = +1(— 1) for g^g^^ > 0(< 0). The positive ApB observed at Tevatron prefers 
to have ^ = — 1. The spin-correlation matrix has 

CO - -^ + ^^' (48) 

= 9-9^' ^''^ 

9(2m + l)2 9A2(2m + l)' ^ ' 

(2m - 1) CQ* 2^sml s^£c0* 
= 9 9A^ + -9A^' ^^^^ 

Diagonalize the C matrix, one has the three eigenvalues as 

, ^\ 2-/32(l-c^.) 2Csce,/3 , (3^s\l + 4.) 

18 ^ A^''' 18 9A2 + 18 A4 " ^^^^ 
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One can easily check that the last eigenvalue is always the largest one, independent of s and A. So, this 
provides us a fixed axis to quantize the spins of top and anti-top quarks. Simplifying the eigenvector 
corresponding to the largest eigenvalue, we have the following best spin quantization axis for spin 
correlation 



(X p + 



ce*{j - 1) 



A2 



k. 



(53) 



with 7 = 1/v^l - /32. For the case that A S> s and the new physics is decoupled, we recover the 
formula for the axis maximizing spin-correlation in the SM [35]. 

Coming back to the axigluon model and neglecting the resonance width, we obtain the best quan- 
tization axis to be 



oc p + 



eg. (7 - 1) 



M2, 



k. 



A. 3 Flavor violating W models 

Considering the following interactions, 

for the interference term, we have 

s {9a + 9v) ilm'^ ((-1 + cf + s + ((1 + cf + W) M^,,) 



lint 



p 



tint 



int 
Co 



int 
C4 



9M2,,(-i + M2„) 
AmsgAgy + c)m'^s - (2 + c)M^,) 

9M^,{-t + M^,) 
2s(3gA9v i^m^ ((-1 + cf +2m)s- ((1 + cf + 2m) M^,) 
9(1 + 2m)M^, {-t + M^,) 

(-1 + c2 + 4m2) s {g\ + gl) 



2 ^ ^"W 



+ 



{g\ + gl) 



rn^s 



2M^,) 



int 
^6 



9(1 + 2m)M^, {-t + M^,) 

(32^3 + 48m4 + Am^ (l + /^^ + 2z'^l3^)) {g\ + g^) {m^s + 2M^,) 
72m2(l + 2m)^M^, {t - M^,) 
2s/3 {gA + 9v) HrrP (c + m + 2m'^) s - {-c + m + 2m2) M^,) 
9(l + 2m)M2„ i-t + M^,) ■ 



(54) 



(55) 

(56) 
(57) 
(58) 

(59) 

(60) 



(61) 
(62) 
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For the new physics squared diagram, 

^ ((5 - 8m^ + c(2 + z/3)) {g^y + g\) + 2 (-1 + Sm^ + 3c(2 + ffjff^) 

8(i-M2„)2 



"i's' (5i + 5^) ' ((-1 + c)2s + 16M2,,) 



4„3 /'„2 

— t— I /t r I II I — t— I I .-1 — t— I I I f I/ / >- t I 

(63) 



ms^gAgv {gl + al) (4^ + M^,) (1(-1 + c)m2s - (1 + c)Ml,) 

^„ew ^ g'/gg^gy (gj + g^) + c)m4(l - c + 2m)g2 _ 2cm3gM^, + (1 + c)(l + c + 2m)M^,) 

^ 2(l + 2m)M^, (-t + M2„)2 

(65) 



new 
Co 



m 



(66) 
(67) 



2g3 (2 ((-1 + c)(l + 3c) + 8m^) g\gl + (-3 + c(2 + c) + Sm^) (g^ + ^4 ^) 

8 {-tMw + M^,) 2 

m's^ {2m^s^g\gl + 2m') s {g\ + Qg\gl + g^) M^,, + (g^ + 5^2 ^2 ^ ^4 ) ^^^4^^ 

2M^, {-t + M^,)' 

WT^^^nin^^^'^ (2™^(1 - c + 2..)^.^gig^+ 

^sM^, (2 (16(1 + c)m^ + 48m^ + (-1 + c)(l + 3c) (-1 + + 8m2 (-1 + 3c2 + /J^)) ^^2 ^2^+ 
(-16(-3 + c)m^ + 80m^ + (-1 + c)(3 + c) (-1 + (3^) + 8^2 (-1 + (l + z') ^')) {g\ + gy)) + 
M^, (2 (5 + 4(7 + 3c)m + 36m2 + /3(3(2 + c)z + 2/3)) g^g^ 

+ (-l + 2c + 4cm-4m(3 + 5m) + (-2 + z2)/32) (^4 ^^4^))^ (68) 
"""'^ ' -2m\l-c + 2m)s'g\gl + 



A{l + 2m)M^,{-t + Ml,)' 
ms (-m(l + 2m) {g\ - g^)^ - c(l + m) {g\ + 6gig^ + g^)) M^,+ 

(l + c + 2m)(g^ + 6gig^ + g^)M^,). (69) 
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